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. Abstract. The paradigm for Anomalous X-ray Pulsars (AXPs) has evolved recently with the dis- 

' covery by INTEGRAL and RXTE of flat, hard X-ray components in three AXPs. These non-thermal 

\ spectral components differ dramatically from the steeper quasi-power-law tails seen in the classic X- 

■ ray band in these sources, and can naturally be attributed to activity in the magnetosphere. Resonant, 
(— I ' magnetic Compton upscattering is a candidate mechanism for generating this new component, since 
Q^, it is very efficient in the strong fields present near AXP surfaces. In this paper, results from an inner 

I ' magnetospheric model for upscattering of surface thermal X-rays in AXPs are presented, using a 

O I kinetic equation formalism and employing a QED magnetic scattering cross section. Characteris- 

tically flat and strongly-polarized emission spectra are produced by non-thermal electrons injected 
in the emission region. Spectral results depend strongly on the observer's orientation and the mag- 
netospheric locale of the scattering, which couple directly to the angular distributions of photons 
sampled. Constraints imposed by the Comptel upper bounds for these AXPs are mentioned. 
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O ; INTRODUCTION 

^' : 

^ . A topical focus of the high energy astrophysics of compact objects over the last decade 

\ has been the so-called magnetars (Duncan & Thompson 1992), constituted by Soft- 

^ \ Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs), whose amassed ob- 

servational properties have indicated that they are isolated neutron stars with ultra- strong 
magnetic fields. The AXPs, are a group of 6-7 pulsating X-ray sources with periods 
\ around 6-12 seconds. They are bright, possessing peak luminosities Lx ~ 10^^ erg s^^ , 

show no sign of any companion, are steadily spinning down and have ages T ^ 10^ 
years (e.g. Vasisht & Gotthelf 1997). Details of the persistent pulsed X-ray emission for 
AXPs are discussed in Tiengo et al. (2002), for XMM observations of IE 1048.1-5937, 
and Juett et al. (2002) and Patel et al. (2003), for the Chandra spectrum of 4U 0142-1-61. 
This emission displays both thermal contributions, which have kT ~ 0.5 — 1 keV and 
so are generally hotter than those in isolated pulsars, and also non-thermal components 
with steep spectra that can be fit by power-laws dn/dE oc E " of index in the range 
5 = 2 - 3.5 (see Perna et al., 2001, for spectral fitting of ASCA data on AXPs). 

The recent detection by the IBIS imager on INTEGRAL, and the PCA and HEXTE 
detectors of the Rossi X-ray Timing Explorer, of hard, non-thermal pulsed tails in three 
AXPs has provided a new twist to the AXP phenomenon. In all of these, the differential 
spectra above 20 keV are extremely flat: IE 1841-045 (Kuiper, Hermsen & Mendez 



2004) has a power-law energy index of s = 0.94 between around 20 keV and 150 keV, 
4U 0142+61 displays an index of * = 0.2 in the 20 keV - 50 keV band, with a steepening 
at higher energies implied by the total DC+pulsed spectrum (Kuiper et al. 2006), and 
RXS J1708-4009 has s = 0.88 between 20 keV and 150 keV (Kuiper et al. 2006). The 
tails are much flatter than the non-thermal spectra in the < 10 keV band, and do not 
continue much beyond the IBIS energy window: there are strongly constraining upper 
bounds from Comptel observations of these sources that necessitate a break somewhere 
in the 150-750 keV band (see Figures 4, 7 and 10 of Kuiper et al. 2006). 

This paper summarizes results from our initial exploration (Baring & Harding 2007) 
of the production of non-thermal X-rays by inverse Compton heating of soft, atmo- 
spheric thermal photons by relativistic electrons, serving as a model for generating the 
hard X-ray tails in AXPs. The electrons are presumed to be accelerated along either open 
or closed field lines with super-Goldreich-Julian densities, perhaps by electrodynamic 
potentials, or large scale currents associated with twists in the magnetic field structure 
(e.g. Thompson & Beloborodov 2005). In the strong fields of the inner magnetospheres 
of AXPs, the inverse Compton scattering is predominantly resonant at the cyclotron fre- 
quency, with an effective cross section well above the classical Thomson value. Hence, 
proximate to the neutron star surface, in regions bathed intensely by the surface soft 
X-rays, this process is extremely efficient for an array of magnetic colatitudes, probably 
dominating other processes such as synchrotron and bremsstrahlung radiation that are 
employed in the models of Thompson & Beloborodov (2005) and Heyl & Hemquist 
(2005). This prospect motivates the investigation of resonant inverse Compton mod- 
els. Here, the general character of emission spectra is presented, using collision integral 
analyses that will set the scene for future explorations using Monte Carlo simulations. 

RESONANT COMPTON UPSCATTERING IN AXPS 

In devising any radiation emission model to describe the non-thermal X-ray tail lumi- 
nosity from AXPs, it is necessary to ascertain the criteria that must be satisfied in order 
to explain the energetics. These were discussed in Baring & Harding (2007), who de- 
termined that for radiative processes that were electromagnetic in origin, i.e. involved 
electrons, the requisite electron densities must be super-Goldreich-Julian. This require- 
ment is rather general in nature, not being constrained to just the resonant Compton 
upscattering scenario that is the focus here, but being coupled to a presence of relativis- 
tic electrons moving along B, with an abundance that can power the intense AXP X-ray 
luminosities, Lx ^ 10^^ erg/sec above 10 keV (Kuiper et al. 2006). The hard X-ray tail 
luminosities are 2-3 orders of magnitude greater than the classical spin-down luminos- 
ity £'sD ~ 87t^B^^ / (SP^c^) due to magnetic dipole radiation torques, where Bp is the 
surface polar field strength, P is the pulsar period and R is the stellar radius. 

Here we briefly recapitulate the energetics analysis of Baring & Harding (2007). 
Let He be the number density of emitting electrons, {je) be their mean Lorentz fac- 
tor, and grad be the radiative efficiency during their traversal of the magnetosphere, ei- 
ther along open or closed field lines. Then Lx ~ ei.ad(7e)'?^eC^(4;rnei?^c) if the emis- 
sion column has a base that is a spherical cap of radius Rc < R. If Rc ~ 10^ cm, 
this yields number densities ~ 3 x 10^^Lx,35/erad(7e) cm^^ for scaled luminosities 



^x,35 = Lx/lO-'^ erg/sec. Comparing eue to the classic Goldreich- Julian (1969) density 
Pgj = V.E /Ak = —D..B/ (iTlc) for force-free, magnetohydrodynamic rotators, one can 
establish the ratio 

en, _ 4,670 LX35P 



(1) 



for AXP pulse periods P in units of seconds, polar magnetic fields Bi^ in units of 
10^^ Gauss, and cap radii R(, in units of 10^ cm. For erad(7e) ^ 10^ and i?6 ~ 1 , the 
requisite density n, is super-Goldreich- Julian, but not dramatically so. No choice of 
radiation mechanism has been invoked in this line of reasoning. Yet observe that electron 
Lorentz factors of ^ 10^ — 10^ and their efficient resonant Compton cooling (i.e. 
£rad ~ 0.01 — 1 ) are readily attained in isolated pulsars with 5 ~ 0.1 (e.g. see Sturner 
1995; Harding & Muslimov 1998; Dyks & Rudak 2000). Such efficiencies should persist 
into the magnetar regime. To substantiate this assertion, observe that in the Thomson 
cooling regime, the resonant Compton cooling rate \%\ ^ (n^cTTc) ZnB^ / (4af7ee^) can 
be obtained from Eq. (17) of Baring & Harding (2007), correcting a missing factor of 
OCf = / {he) in the denominator. The lengthscale Ares — c/\Ye\ for resonant cooling can 
then be estimated from the Planck spectrum photon number density Ug ~ S^'^ j'k'^ at the 
surface, where the scaled temperature ^ — kT /{mgC^) sets ~ 3^ , and Ac = h/nigC 
is the electron Compton wavelength over lit . The result is Arcs ~ (Ac/cUt) Jaj^^B^) , 
lengthens considerably at high altitudes due to dilution of the soft photon density. For 
5~0.1, 10^ and ST ~ 3 X 10 ^ , Ares is much less than the stellar radius, a result 
that is clearly not modified by relativistic quantum corrections when 5^1. 

By necessity, the locales where Compton interactions sample the cyclotron resonance 
are confined to the lower altitudes in an AXP magnetosphere, where the field is suffi- 
ciently high. This is controlled primarily by the scattering kinematics, which dictates a 
coupling between the energies CymgC^ and 7emeC^ of colliding X-ray photons and elec- 
trons, respectively, and the local angle Qy of the interacting photon to the magnetic field 
lines. The cyclotron fundamental is sampled when 

7eey(l -cos0y) fa B , for7e > 1 . (2) 

For X-ray photons emanating from a single point on the stellar surface. Baring & 
Harding (2007) computed the zones of influence of the resonant Compton process for 
dipole field geometry in a flat spacetime. Assuming that the X-rays propagate with 
no azimuthal component to their momenta, the resonance criterion is satisfied on a 
surface that is azimuthally symmetric about the magnetic field axis, for slow rotators. 
For outward-going electrons, the locus of the projection of this surface onto a plane 
intersecting the magnetic axis was found to be 

3 _ m a/1 + 3cos^0 ^ _ Bp 

^ l-COS% ' Ije&y ' 

where x = r/R is the altitude scaled in units of the neutron star radius R, and 6 is the 
magnetic colatitude of the point of scattering. Bp is the surface polar field strength in 

units of 5cr = 4.413 x lO^-' Gauss. Geometry determines the function, dy = 0y{6e, 0) , 
which is given in Eq. (5) of Baring & Harding (2007), where dg is the colatitude of 
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FIGURE 1. Contours in a section of a pulsar magnetosphere that depict cross sections of the surfaces of 
last resonant scattering by ultra-relativistic electons, i.e. the maximal extent of the Compton resonasphere. 
The heavyweight blue contours are computed for different values of the resonance parameter defined 
in Eq. (3), and at extremely high altitudes asymptotically approach the magnetic axis (central vertical 
line). The filled red circles denote the neutron star, whose radius R establishes the linear spatial scale for 
the Figure. The cases illustrated are for photons emanating from the polar axis (i.e. 9e —0° , left panel) 
and from the equator (i.e. Og = 90° right panel), denoted by the small black dots, for which the neutron 
star shadow regions are demarcated by the light dotted boundaries, and only the surfaces (azimuthally- 
symmetric about the magnetic axis) are accessible to resonant Compton interactions. 



the surface emission point. Here ^ is the key parameter that scales the akitude of the 
locale of resonant interaction, and typically falls in the range 1 — 10^ for magnetars 
when 7e ~ 10^ — 10"^ . For the broadly representative cases of soft photons emitted from 
the surface pole (Oe = 0) and magnetic equator (6e = 90°), the surfaces of resonant 
scattering for different ^ are illustrated in Fig. 1. Shadows of the emission points are 
also indicated to mark propagation exclusion zones for the chosen emission colatitudes. 

It is evident from the Figure that the altitude of resonance is much lower for equatorial 
emission cases (comparing left and right panels), and also in the equatorial regions 
when compared with polar locales (within each panel). At small colatitudes above 
the magnetic pole, Oy is necessarily small, pushing the resonant surface to very high 
altitudes where the field is much lower. In equatorial interaction locales, which are 
preferentially sampled for quasi-equatorial emission colatitudes, the photons tend to 
travel more across field lines in the observer's frame, and so access the resonance in 
regions of higher field strength, thereby reducing the altitudes where the resonance is 
sampled. These are manifestations of the correlation between Gy and B (for fixed Ey 
and Ye) evinced in Eq. (2). For cases < 0^ < 90° (not depicted), the contours are 
morphologically similar, though they incur significant deviations from those in Fig. 1. 
Clearly, by sampling different emission colatitudes Oe these surfaces are smeared out 



into annular volumes. Observe that introducing an azimuthal component to the photon 
momentum tends to increase propagation across the field, i.e. raising Oy, so that the 
resonance is accessed at lower altitudes and higher field locales. Hence, loci like those 
depicted in Fig. 1 actually represent the outermost extent of resonant interaction, and 
so are surfaces of last resonant scattering, i.e. the outer boundaries to the Compton 
resonasphere. It is evident that, for the majority of closed field lines for long period 
AXPs, this resonasphere is confined to within a few stellar radii of the surface. 



Resonant Compton Upscattering Spectra 

Collision integral calculations for upscattering spectra from resonant Compton inter- 
actions are routinely obtained for uniform magnetic fields. In the AXP problem, the 
scalelength Arcs for the interaction is often much shorter than the scale R of the field 
divergence/gradient, so such uniform B computations are reasonably informative. Let 
ny be the number density of photons resulting from the resonant upscattering process. 
For inverse Compton scattering, an expression for the spectrum of photon production 
dny/{dtdefdilf) , differential in the photon's post-scattering laboratory frame quanti- 
ties Ef and Hf — cos©/, was presented in Eqs. (A7)-(A9) of Ho and Epstein (1989), 
valid for general scattering scenarios. The dimensionless pre- and post-scattering photon 
energies (i.e. scaled by nigC^ ) in the observer's frame (OF) are and £/ , respectively, 
and the corresponding angles of these photons with respect to the OF electron velocity 

vector — j8e (i.e. field direction) are 0, and 0/ , respectively. The result for the spectrum 
can be integrated over fif and then written as (detailed in Baring & Harding, 2007) 



dny _ ngns c 



dtdSf jU+ — jU_ 



Here, the notation /i,- = cos0, and jif = cos0/ is used for compactness, n^ and «^ 
are the number densities of relativistic electrons and soft photons, respectively. For 
simplicity, the incident photons are assumed to be monoenergetic and to possess a 
uniform distribution of angle cosines jU/ in some range /i_ < jU,- < /i+ . Note also 
that since ^ I , the incident photon angle with respect to B in the electron rest 
frame (ERF) is 0/ fti 0, while the ERF angle 6f of the scattered photon relative to 
the field can take on a range of values. The function tt)'(tt),, Of), which appears in 
the 5 function in Eq. (4), encapsulates the electron rest frame scattering kinematics, 
as detailed in Baring & Harding (2007). The differential cross section, do /d{cos6f) , 
appearing in Eq. (4) is evaluated in the ERF, and is taken from Eq. (23) of Gonthier et 
al. (2000); it incorporates relativistic QED physics that is applicable for arbitrary field 
strengths. Specialized to the case of scatterings that leave the electron in the ground 
state, the zeroth Landau level that it originates from, its specific form is summarized 
in Baring & Harding (2007). More general results for the fully relativistic, quantum 
cross section for resonant Compton scattering can be found in Herold (1979), Daugherty 
& Harding (1986), and Bussard, Alexander & Meszaros (1986). These extend earlier 
non-relativistic quantum mechanical formulations such as in Canuto, Lodenquai & 
Ruderman (1971), and Blandford & Scharlemann (1976). 



Resonant Compton Upscattering Spectra 
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FIGURE 2. Resonant Compton upscattering spectra (scaled, unpolarized) such as might be sampled 
in the magnetosphere of an AXP, for different relativistic electron Lorentz factors /g , as labelled. The 
emergent photon energy Ef is scaled in terms of m^c^ . The chosen magnetic field strengths of Z? = 3Bcr 
(heavyweight, blue) and B = 0.3Bcr (lighter weight, red) correspond to different altitudes and perhaps 
colatitudes. Results are depicted for seed photons of energy = 0.003 (marked by the green vertical line), 
typical of thermal X-rays emanating from AXP surfaces; downscattering resonant emission at Ef < Ej 
was not exhibited. Specific emergent angles of the emission in the observer's frame, with respect to the 
magnetic field direction, are indicated by the filled magenta symbols for four of the spectra, with triangles 
denoting &f = 5° , squares corresponding to ©/ = 1° , and circles representing 0/ = 0.2° . 



The relativistic Compton cross section o is strongly peaked at the cyclotron funda- 
mental (see Fig. 2 of Gonthier et al. 2000) due to the appearance of resonant denomi- 
nators in the S-matrices. These generate a Lorentz profile factor l/[((0/ — 5)^ + (r/2)^] 

in do / d {cos, 9 f) , where T = T^y^Ei and Ei = ^ I + 2B + aif is the energy of the in- 
termediate electron. Here F^yc <^ 5 is the dimensionless cyclotron decay rate from the 
first Landau level (with electron momentum = CO,- parallel to B), signifying that the 
intermediate electron states become effectively real at the cyclotron resonance, and pos- 
sess a finite decay lifetime. A form for the decay rate Fcyc , can be found in Eqs. (13) or 
(23) of Baring, Gonthier & Harding (2005; see also Fatal 1986; Harding & Fai 2006). 
For 5^1, Fcyc ~ 4af5^/3, while for 5^1, quantum and recoil effects generate 
Fcyc ~ (0Cf/e)A/5/2. Note that Baring & Harding (2007) adopted the ansatz F Fcyc 
appropriate for resonant Thomson scattering; here this is updated to incorporate rel- 
ativistic corrections to the resonance width following Harding & Daugherty (1991), 
amounting to a spectral normalization correction by a factor of 1 / \/l+2B for all ef . 

Representative spectral forms are depicted in Fig. 2, for the situation where emergent 
polarizations are not observed (see Baring & Harding 2007, for polarization characteris- 
tics). Because of the narrowness of the resonance, non-resonant scattering contributions 



were omitted when generating the curves; these contributions produce steep wings to 
the spectra at the uppermost and lowermost energies, and a slight bolstering of the flat 
portion. Generally they contribute significantly only when access to the resonance is 
kinematically forbidden, i.e. outside the Compton resonaspheres illustrated in Fig. 1. 
The resonant restriction kinematically limits the emergent photon energies £/ to 

Ye{l-Pe)B <ef< ^^^^^^ , (5) 

a range that generally extends below the thermal photon seed energy . For the 5 = 0.3 
case in Fig. 2, a quasi-Thomson regime, the spectra are characteristically flat (e.g. see 
Dermer 1990; Baring 1994; Liu et al. 2006) for most £/, indicative of the kinematic 
sampling of the resonance in the integrations over soft photon angles 0, . Since F <C 5 in 
general, the normalization of this flat portion scales as driy/ {dtdEf) oc l/F. Only at the 
highest energies does the spectrum begin to deviate from flat (i.e. horizontal) behavior, 
and this domain corresponds to significant scattering angles in the ERF, i.e. cosines 
1 — cos Bf not much less than unity. Then the mathematical form of the differential 
cross section becomes influential in determining the spectral shape, as discussed in 
Baring & Harding (2007). In AXPs, the B = 0.3 case best represents higher altitude 
locales for the resonasphere, such as at smaller colatitudes near the polar axis. Fig. 2 
also exhibits spectra for B = 3 , a case more typical of equatorial resonance locales. Then 
the flat spectrum still appears at energies < £/ <C 7e(l + j8e)5/(l + IB) , when again 
cos dfTtil. Yet the 5 = 3 curves in the Figure display more prominent reductions at the 
uppermost energies £/ ~ /^(l + j8e)5/(l +25) due to the sampling of 1/2 < coy <^ co, 
values in the ERF that correspond to strong electron recoil effects. 

Photons emitted in these uppermost energies have 1 — cos 0/ ~ 1 in the ERF and are 
highly beamed along the field in the observer's frame. This is an important property 
that is highlighted via the filled symbols in Fig. 2. The intense beaming of radiation 
along B in the OF, and the profound correlation of the angle of emission @f with 
the emergent photon energy £/, are both consequences of scattering kinematics in 
the resonance. In the resonant case here, the 7e ^ 1 regime dictates that most of 
the emission is coUimated to within 5° of the field direction, and rapidly becomes 
beamed to within 0.2° as the final photon energy increases towards its maximum. This 
kinematic characteristic guarantees that spectral formation in Compton upscattering 
models is extremely sensitive to the observer's viewing orientation in relation to the 
magneto spheric geometry. This suggests powerful geometrical probes of AXP emission 
regions if pulse-phase spectroscopy is achievable in future generation observatories. 

CONCLUSION 

This paper has summarized some essentials of the resonant Compton upscattering model 
for hard X-ray tail emission in AXPs, as developed in Baring & Harding (2007). The 
spectra exhibited in Fig. 2 are considerably flatter than the hard X-ray tails ( ~ e^^" for 
a ~ 0.2 — 1 ) seen in the AXPs, and extend to GL AST-band energies much higher than 
can be permitted (i.e. around 750 keV) by the Comptel upper bounds to these sources, 
unless emergent angles 0/ to B exceed around 1°. They represent a preliminary in- 



dication of how flat the resonant scattering process can render the spectrum, which can 
readily be steepened by spatial distribution of electron injection, significant and unavoid- 
able cooling, and also non-resonant contributions. What an observer detects will depend 
critically on his/her viewing perspective and the magnetospheric locale of the scattering. 
The one-to-one kinematic coupling between ef and jUy implies that the highest energy 
photons are beamed strongly along the local field direction. This may or may not be 
sampled by an instantaneous observation at a given rotational phase. Realistically, for 
many pulse phases, angles corresponding to /i/ < 1 will be predominant, lowering the 
value of Ef . How low is presently unclear, and remains to be explored via a model with 
full magnetospheric geometry, the obvious next development in this research program. 
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